Investigation of the ^^Na nucleus via resonance elastic scattering 
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Abstract 

The structure of the unbound proton-rich isotope ^^Na was studied in resonance elastic scattering 
of a radioactive ^^Ne beam on a proton target using the thick-target inverse-kinematics method. 
The experiment covered from 0.5 meV to 2.7 meV in c.m.s. Only one state of ^^Na (the second 
excited state) was observed. A combined i?-matrix and potential model analysis was performed. 
The spin and parity assignment of this second excited state was confirmed to be 1/2^. We show 
that the position of the 1/2"'" state significantly affects the reaction rate through that state but the 
total reaction rate remains unchanged since the ^^Ne(2p,7) proceeds mostly via the ground and 
first excited states in ^^Na at stellar temperatures. 
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I. INTRODUCTION 



The structure of light exotic isotopes is currently one of the major topics in nuclear 
physics. Properties of light neutron-deficient isotopes are not as well known as the properties 
of their neutron-rich mirrors. The Coulomb interaction leads to proton unbound states in 
neutron-deficient nuclei even though the corresponding states are bound in the mirror nuclei. 
Decay properties of these nuclei yield additional information about their structure as well 
as the structure of states in the mirror nuclei. One such example is ^^Na, a proton drip-line 
isotope which is unbound with respect to proton emission by only 321 keV. It was first 
observed by Cerny et al. ^ and until very recently only its mass and the excitation energy 
of its first excited state were experimentally known. The lack of experimental information 
on the level structure of ^^Na is due to technical difficulties related to the population of 
the excited states of this rather exotic nucleus with three fewer neutrons than the stable 
sodium isotope ^^Na. The second excited state in ^^Na was recently observed in a ^^Ne+p 
elastic scattering experiment [3] in which a narrow excitation energy range (0.4 - 1.1 MeV) 
was covered. The original goal of this work was to extend our knowledge of the structure 
of ^^Na to higher excitation energies (up to 2.7 meV) using the resonance elastic scattering 
of a radioactive beam of ^^Ne on protons (the upper limit in excitation energy is implied by 
the maximum possible energy of ^^Ne beam that can be produced at ISNAP P]). During 
the review stage of this manuscript new study of the ^^Na nucleus was published from the 
study at SPIRAL (GANIL) by Oliveira et al. Q]. In this work, the spectrum of ^^Na was 
measured up to excitation energy of 5.7 MeV at one angle (180°) using the Thick Target 
Inverse Kinematics (TTIK) technique j^. Two additional resonances at higher excitation 
energies were found in In addition, two broad peaks at Ecm ~ 2.4 and 3.1 MeV were 
observed. It was pointed out by the authors in 4J that these peaks can originate from 
process other than elastic scattering. The first peak is in the energy range of our experiment 
and should be populated if the origin of this peak is in fact related to elastic scattering. 
Therefore, it is of interest to compare our data, obtained at a lower beam energy and 
different angles with the data from 4]. The structure of low lying resonances in ^^Na also 
has some astrophysical interest. It was shown in that at certain stellar temperature 
and density conditions two-proton capture reaction on ^^Ne can play an important role in 
bridging the waiting point nucleus ^^Ne and provide continuous flow between the CNO and 
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the Fe-Ni mass region in the rp -process. The ^^Ne(2p,7) reaction rate was recalculated 
taking into account new information on the excitation energy of the second excited state. 

In the following sections, we present the results of resonance elastic scattering of ^^Ne on 
protons, discuss the structure of the ^^Na nucleus, and finally discuss the influence of a 1/2"'" 
state on two-proton capture reaction rate on the ^^Ne nucleus. 

II. EXPERIMENT 

The experiment was carried out at the TwinSol radioactive nuclear beam facility Q of 
ISNAP (Institute of Structure and Nuclear Astrophysics) at the University of Notre Dame. 
A beam of ^^Ne was produced via the ^He(^^0,^^Ne)n reaction. The experimental setup is 
shown in Fig. ^ A primary beam of ^^O with an average intensity of 230 electrical nA and 
an energy of 80 MeV was incident on a 2.5 cm long gas cell containing '^He gas at a pressure 
of 1 atm. A Faraday cup placed after the gas cell was used to stop the primary beam and to 
measure its intensity. Two superconductive solenoids acted as thick lenses to separate the 
^^Ne beam from other reaction products. The magnetic rigidity of the first superconductive 
solenoid was set to focus ^^Ne"'"^ ions with an energy of 57 MeV onto a 10 fig/ cm^ C stripping 
foil as shown in Fig. ^ while the second solenoid focused fully-stripped ^'^Ne"'"^° ions onto 
the secondary target. The different rigidities of the two solenoids allowed for a much better 
separation of ^^Ne from other reaction products. The ^*Ne ions were focused onto a 1 cm 
(diameter) spot at the secondary target position. Under these conditions, a ^^Ne beam was 
obtained with an intensity of about 10^ pps, an energy of 56.3 MeV, and an energy spread 
of 3.2 MeV full width at half maximum (FWHM). The desired ^^Ne ions made up only 
2% of the secondary beam; the main contamination was ^^O (see Fig. Polyethylene 
(CH2) foil was used as the secondary target. Its thickness of (5.52 mg/cm^) was chosen to 
stop the ^^Ne ions completely. The Thick- Target Inverse-Kinematics (TTIK) technique 
was used to obtain an excitation function for resonance elastic scattering. This method is 
based on the large differences of specific energy losses for heavy ions (^^Ne) and protons. 
Protons backscattered by ^'^Ne penetrate the target and lose only a small fraction of their 
energy. The center-of-momentum energy at the moment of the interaction can therefore be 
recovered from the measured proton energy. Three Si detectors placed at angles of 7.5°, 
22.5° and 37.5° with respect to the beam axis (corresponding to backward angles for elastic 
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FIG. 1: The experimental setup with beam trajectories through TwinSol. The "lohipop" reduces 
contamination of the beam by intercepting ions that focus at a different location relative to the 
""^^Ne beam. 

scattering of ^^Ne on protons) were used to detect recoil protons. The primary ^^O beam 
was pulsed with a repetition rate of 5 MHz and a pulse width of approximately 2 ns. The 
5.6 m distance from the primary target to the array of silicon detectors was adequate to 
obtain a clean separation of the protons associated with the interaction of ^^Ne ions in the 
target from protons associated with elastic scattering of beam contaminants and reactions 
of other beam isotopes in the target (see Fig. Spectra of protons (in the laboratory 
frame) measured in the three detectors are shown in Fig. El Measurements with a thick 
natural carbon target were performed to subtract the background associated with carbon 
in the polyethelyne target. The thickness of the carbon target was matched to that of the 
polyethylene target in terms of proton energy losses. The shaded areas in Fig. |21 show the 
contributions from reactions of ^'^Ne on the carbon target. 

Proton spectra taken at three different angles, converted to the cm. system, are shown 
in Fig. ^ In the conversion of raw spectra to the absolute cross section we used the mea- 
sured ratio of ^^Ne to the intensity of the primary beam ^^O. The kinematic relations and 
energy losses of protons in the target were taken into account in the transformation to the 
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FIG. 2: Composition of the radioactive beam measured by a Si detector placed at the secondary 
target position. The energy of the ions is shown versus the time of flight (ToF). The "Start" signal 
for the ToF measurements was produced by the Si detector, and the "Stop" signal was derived 
from the radiofrequency (RF) signal of the buncher. 

cm. excitation functions. Monte Carlo simulations of the experimental energy resolution 
were performed for several excitation energies and angles. The typical experimental reso- 
lution at 165° was found to be about 30 keV at cm. energies above 1 MeV. The precision 
of the excitation energies is about 20 keV, defined mainly by the uncertainties in the spe- 
cific energy losses of ^^Ne ions in the CH2 target. The absolute cross section values are 
obtained with uncertainties of 15%. Only one resonance was observed in the ^^Ne+p elastic- 
scattering excitation function in the measured energy range up to 2.7 MeV. This resonance 
was previously reported in [2] and j^. We confirm the 1/2"*" spin-parity assignment for 
the resonance. The interference pattern of the / = resonance phase-shift with Coulomb 
phase-shift and the angular distribution (three different angles were measured) provides for 
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unique spin-parity assignment 1/2"'" for the observed resonance. The excitation energy and 
width of this state were obtained in an /^-Matrix fit: Eex= 0.74 ± 0.03 MeV and F = 130 
± 50 keV. Statistical uncertainties dominate the reported errors. The parameters of the 
resonance are in agreement with the results obtained for this state in [2]. Continuum Shell 
Model (CSM) calculations [7J, performed for a chain of oxygen isotopes, predict that the 
1/2+ second excited state in ^^O (the mirror of the observed state in ^^Na) has an almost 
pure (~ 70%) single particle nature. We need a proper determination of the single parti- 
cle spectroscopic factor for this resonance to test the predictions of the CSM. Often, the 
Wigner limit 7^ = is used to estimate the single particle width of the resonance and 
to extract a spectroscopic factor. However, this approach gives only an order-of-magnitude 
estimate. When applied to the observed 1/2"*" state, a misleading result is obtained. In 
)articular, the single particle spectroscopic strength extracted for the state is only 25% 
3]. In this work, we use the potential- model approach instead. In this case, the properties 
of the ^^Ne+p system are described by a spherical Woods-Saxon potential. If the param- 
eters of the potential are known, then one can predict the properties of the single particle 
states. The width of the states calculated with the proper potential can then be compared 
with the experimental values and the spectroscopic factor can be extracted as the ratio of 
experimental width to the calculated value (as long as only one decay channel is open, as 
is the case in the energy region of interest). The same approach was adopted in P] but no 
details on the potential parameters are given. The initial parameters of the Woods-Saxon 
potential were found by fitting the excitation energies of the lowest 5/2"*", 1/2"^, and 3/2"*" 
levels in ^^O and ^^F simultaneously It is known that these states are pure single particle 
in nature. (The Spectroscopic Factors (SF) are 0.9, 0.9 and 1.25, respectively, ^j). A single 
set of Woods-Saxon potential parameters generates correct energies for the ground and the 
first two excited states in both ^'^O and ^'^F to within 100 keV. The optimum parameters are 
shown in the table HI The width of the 3/2+ state (the only unbound state among the three 
in ^^0/^'^F) calculated with the potential given in the tablelJis 90 keV for ^^O and 1.14 MeV 
for ^''F, Q]. The experimental values are 95 keV and 1.5 MeV, respectively. The ratio of the 
experimental to the calculated width gives the spectroscopic factor, which is greater than 
one for both ^''O and ^''F. To avoid this, we could increase the diffuseness of the potential 
( as it was suggested in jsi]) to get SF=1.0 for 3/2+ instead of 1.25. But since we want to 
keep all potential parameters the same for ^^0/^^F and ^^0/^^Na (except Vq), we decided 
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TABLE I: Parameter of the Woods-Saxon potential 



Parameter 


170 /i^F 


i90/i9Na 


^0 


-55.77 


-50.12 


Vosl 


6.40 


6.40 


ro 


1.21 


1.21 


rosi 


1.21 


1.21 


rcoui 


1.21 


1.21 


a 


0.65 


0.65 


a-sl 


0.65 


0.65 



to use standard values. A potential with the same Woods-Saxon parameters (except for Vq 
parameter) was used to calculate the widths and excitation energies of the single particle 
states in ^^O. The potential was adjusted to fit the position of the first 1/2^ state in ^^O 
which has pure single particle structure (see Table H}. The excitation energy of the 
other pure single particle state (3/2+at 5.45 MeV) was reproduced by this potential with 
an accuracy of 40 keV. The potential model predicts the width of this latter state to be 330 
keV, which is consistent with the experimental value of 320 ± 25 keV 10]. The ground state 
of ^^O appears to be underbound by 500 keV, but it is known from shell model calculations 
0] and also from experiment that the ^^O ground state has significant admixtures from 
configurations other than simple ((^5/2)^ for the neutrons in excess of A'^ = 8. Thus, our 
potential model is not adequate for a description of the ground state of ^^O. The potential 
obtained in the procedure described above can now be used to calculate the single-particle 
width of the 1/2+ state in ^^Na. It was found that the predicted width of this state is 134 
keV, in comparison with the experimental value of 130 ± 50 keV gives a clear indication 
that it is a pure single particle state (SF = 0.95 ± 0.3). The spectroscopic factor defined in 
this way is in very good agreement with the result of CSM which gives 70% SF and a width 
for the 1/2+ resonance in ^^Na of about 94 keV (obtained as follows: the single particle 
width from potential model was multiplied by SF given by CSM). The excitation function 
for ^^Ne+p elastic scattering can be calculated using the potential with the parameters given 
in the table HI The dashed curve in Fig. |3] shows the result of such a calculation and it is 
almost identical to the i?-matrix fit, shown with the solid curve. 
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FIG. 3: The spectra of protons from the polyethylene target observed in Si detectors placed at 7.5°, 
22.5° and 37.5° in the laboratory frame. The shaded areas represent the proton spectra measured 
with a natural carbon target. The spectra were gated only on the ToF signal of the -'^^Ne-'^'^^ beam 
(see Fig. HJ. 

The same procedure can be used to obtain the SF of the 1/2"*" state from the experimental 
data of [2| by Angulo et al., where the low-energy part {E^rn i 1-4 MeV) of the ^^Ne+p elastic 
scattering excitation function was measured. The width of the 1/2+ state was found to be 
r = 100 ± 8 keV 0]. (The error bars include the uncertainty in the arbitrary choice of 
the channel-radius parameter in the i?-matrix fit). The fit to the data from [2] and our 
potential model gives a spectroscopic factor of 0.75 ± 0.06. The error reflects the statistical 
uncertainty of the data from j^. The state observed at 2.37 MeV in ^^O does not manifest 
itself in the measured excitation function. Nevertheless, conclusions can be drawn from the 
fact that we did not see any indications of the presence of this state in our data. Specifically, 
this state is either very narrow and/or only weakly populated in the elastic channel. The 
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FIG. 4: The excitation functions for Ne+p elastic scattering. The sohd curves correspond to 
i?-Matrix calculations, the dashed curve(only top figure) to potential-model calculations, and the 
dashed-dotted curve to Coulomb scattering. 
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FIG. 5: The reaction rate for two-proton capture ^^Ne(2p,7). The calculations were performed in 
noneresonant and resonant analytical approach used in [6[. The thin, thick and the slashed lines 
overlaps among each other which show reaction rate for two proton capture on ^^Ne which includes 
the contribution from ground state (5/2"'" in -"^^Na ) and the first two excited states (3/2"'", 1/2"'" 
in ^^Na) except slashed line in which excluded contribution from 1/2"*" state in ^^Na. Dashed and 
dot-dashed lines present contribution from the 1/2+ level to the two-proton capture with old and 
new position of 1/2+ levels, respectively. 

CSM calculations indicate that the third excited state of ^^O /^^Na, is built predominantly on 
the first and the second excited states of ^^0/^^Ne, and has a spin-parity of 9/2"*". The width 
of this state in ^^Na, according to the CSM prediction [3] is 1.1 eV. However, this model 
does not take into account possible admixtures from (79/2 shell (only sd shell is included). 
The single-particle width of this state calculated with the potential given in Table |l] is 0.5 
keV. So even if we assume the 10% admixture of the ggj2 shell to the wave function of this 
state the width would become ~ 5 eV. Much too narrow to be measured in the present work. 
Hence, the absence of the state around 2.4 MeV (excitation energy) in the observed ^^Ne+p 
spectrum lends indirect support to the prediction of the CSM regarding the nature of this 
state. 
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In 0, two peaks with cross sections of ~ 300 mb/sr were observed at cm energies of 2.4 
and 3.1 MeV in the excitation function of ^^Ne+p elastic scattering. As shown in Figure 
m no peak is present at 2.4 MeV in our spectrum. It is a clear experimental proof that 
peaks observed in ^ are not related to the "Na states populated in elastic scattering. 
Authors of ^ argue that these peaks can be related to the sequential 2p-decay of the highly 
excited states of ^^Na. Absence of these peaks in our spectrum lends indirect support to this 
hypothesis because due to lower initial ^^Ne beam energy in our experiment these highly 
excited states would not be populated and thus the peaks at 2.4 MeV and 3.1 MeV cannot 
be observed. 

Potential model with potential parameters given in Table |1] reproduces excitation energies 
(with 100 keV accuracy) and widths of all states of single particle nature in ^'^F, ^^O and 
^^O. It also gives a very good description of the ^^Ne+p excitation function up to 2.7 meV 
(cm.) (dashed curve in Fig. |3^). Using this model one can make a prediction for the 
excitation energy and width of the 3/2"'" (ld3/2)^ single particle state in ^^Na. According 
to our calculations it should have an excitation energy of 5.2 meV and width of ~ 2 meV. 
This state is the analog of the known 5.45 meV 3/2"*" single particle state in ^^O P]. It is 
also a feature of both CSM calculations ^ and shell-model calculation using code Oxbash 
jll | with the WBT 12] interactions where it appears at 5.53 meV. Thus this state should 
clearly be in accessible range of where the excitation function was measured up to 5.7 
meV. However, we saw no indication of this state. It would be of great interest to locate 
this state in ^^Na in future measurements or understand the reason for its disappearance. 

It was shown in that for high temperature and density conditions, two-proton capture 
on ^^Ne will compete with {a,p) and decay rate. Such stellar condition would allow a 
fast depletion of ^'^Ne and a fast leakage out of the hot CNO cycle toward higher masses. 
In joj, the second excited state 1/2-1- in the ^^Na nucleus, was not known. The position of 
the 1/2-1- state was estimated from Thomas-Ehrman shift relative to the mirror l/2-|- state 
in ^^O and the value Eg^ = 0.353 meV was used. It is more than two times lower than 
the experimental value obtained by us and as a result some changes in reaction rates niay 
occurs. We recalculated the ^^Ne(2p,7) reaction rate for the same stellar conditions as in pj. 
For the case discussed here, the approach proposed in [l^ can be used. Using the obtained 
energy of the l/2-|- state, the ^^Ne(2p,7) reaction rate was reduced by a factor of 100 due 
to the new location of the 1/2+ resonance (Fig. 13 the ordinate represents total two proton 
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reaction rate on ^^Ne). However, the total reaction rate is almost unaffected since at stellar 
temperatures, the reaction rate is dominated mostly by g.s. and first excited state. 

III. CONCLUSION 

The structure of the proton-rich particle unstable nucleus ^^Na was investigated in the 
excitation energy region from 0.5 meV to 2.7 meV by means of resonance elastic scattering 
of a radioactive beam of ^^Ne on protons. Only an excitation energy of 0.74 

meV, was observed. Its spin and parity were confirmed to be 1/2"'". A potential model 
approach was used to obtain the single particle spectroscopic factor of this resonance. We 
show that the state has a reasonably pure single-particle structure (2sl/2)^, in agreement 
with shell- model calculations The absence of another low-lying state (observed at an 
excitation energy of 2.37 meV in the mirror nucleus ^^O) in the ^'^Ne-l-p excitation function 
is indirect confirmation of the shell- model prediction that this is a 9/2'^ state built on an 
excited states of ^'^Ne. As such, it would be very narrow and only weakly populated in elastic 
scattering. 

We did not observe a broad peak at 2.4 meV, found in jj], which indicates that this 
peak is unrelated to elastic scattering and may be the result of sequential 2p decay from the 
higher lying excited states of ^^Na (as argued in The potential model with potential 
parameters obtained in this work was used to predict the excitation energy and width of a 
3/2"'" single particle state in ^^Na. This state should be well in the range of the experimental 
data of work ^, but there is no indication of this state. Additional measurements at higher 
excitation energies (and different angles) are needed to locate this resonance or understand 
the reason for the disappearance of this state. 

The ^^Ne(2p,7) astrophysical reaction rate was recalculated using the new experimental 
information on the excitation energy of the 1/2"'" state. The (2p,7) reaction rate due to this 
state is actually 100 times lower than it was suggested in previous work 6], however, due 
to the fact that the ^^Ne(2p,7) reaction at stellar temperatures is dominated by the ground 
and the first excited state this finding did not make significant impact on the total reaction 
rate. 
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